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EQUIVALENCE OF ENSEMBLES FOR LARGE 
VEHICLE-SHARING MODELS 

CHRISTINE FRICKER AND DANIELLE TIBI 


Abstract. For a class of large closed Jackson networks submitted to capac¬ 
ity constraints, asymptotic independence of the nodes in normal traffic phase 
is proved at stationarity under mild assumptions, using a Local Limit Theo¬ 
rem. The limiting distributions of the queues are explicit. In the Statistical 
Mechanics terminology, the equivalence of ensembles - canonical and grand 
canonical - is proved for specific marginals. The framework includes the case 
of networks with two types of nodes: single server/finite capacity nodes and 
infinite servers/infinite capacity nodes, that can be taken as basic models for 
bike-sharing systems. The effect of local saturation is modeled by generalized 
blocking and rerouting procedures, under which the stationary state is proved 
to have product-form. The grand canonical approximation can then be used 
for adjusting the total number of bikes and the capacities of the stations to 
the expected demand. 


1. Introduction 

Many cities are now equipped with bike-sharing systems. Those have met a great 
success among the public, but they turn out to face some difficulty for adapting to 
the demand of users. One major inherent problem is due to the spatial heterogeneity 
of the demand. It results in the common observation that, during the day, some 
parking stations stay most of the time empty, while others keep many vehicles 
unused. Some sites may indeed be more popular for picking a vehicle than for 
returning it. This induces two possible drawbacks for users, namely, finding no 
vehicle for rent at their starting point, or no parking room at the end of their 
ride. Indeed, since parked vehicles are locked, each station has only a given number 
of docking places. This capacity limitation may be most problematic, making 
customers loose time in an additional run (and then possibly a walk back to their 
former destination) just for parking their vehicle. On the contrary, finding no vehicle 
at some station is more simply resolved - using alternative means of transport. A 
major concern, when designing such a network, is thus to find the best acceptable 
tradeoff, relating to the total number of offered vehicles, as function of the different 
capacities of the stations. A low supply of vehicles produces empty stations, while 
a large supply causes saturation. 

The purpose of the present paper is to analyze large, inhomogeneous networks 
with finite capacity stations. The performance of such models is estimated through 
explicit asymptotics of their stationary state. Performance indicators such as the 
total rate of failure can then be optimized by adapting the flexible parameters (the 
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total number of vehicles and, possibly, the sizes of the stations). Note, yet, that 
this may not be sufficient to reach a good quality of service. Ultimately, different 
devices should be implemented in order to limit service failure, or disparity between 
stations: For example, setting up a vehicle reservation service, or supplying empty 
stations with vehicles (moved from saturated stations), or else displaying online 
updates of the current state of the network (so that one knows where parking room 
or vehicles are available). In this regard, besides quantifying the weaknesses of the 
simplest primitive system, we introduce state-dependent routing procedures that 
may open the way to such improvements. 

Time inhomogeneity is another concern: The behavior of users should vary in 
time, more or less obeying a 24 hours cycle. This issue will not be addressed in this 
paper, where all processes will be assumed to have constant traffic parameters and 
to have reached their stationary regime. In this respect, a major issue would be to 
estimate the relaxation time - or time to reach equilibrium - of such processes, and 
compare it with the duration of the constant parameter phase of interest. 

It has been commonly observed that vehicle-sharing systems can be modeled 
by closed Jackson networks. In this description, the vehicles play the role of the 
customers - which number is fixed - while the users act as successive servers at 
the parking stations. The associated network has two types of nodes: one server 
nodes, that describe the stations, and infinite server nodes, representing the differ¬ 
ent routes between the stations. In [6] a detailed analysis of large, infinite capacity, 
closed Jackson networks at equilibrium is proposed. Application to vehicle sharing 
systems - at their early experimental stage in the 90’s - is briefly considered. The 
asymptotics of a model mentioned in [6] is examined in [101 as the number of nodes 
(stations and routes) is fixed, while the number of customers (vehicles) grows to 
infinity. Both papers crucially rely on the explicit product-form of the stationary 
distribution, which is standard in the infinite capacity case. 

In 0, assuming complete homogeneity of the traffic and equal finite capacities of 
the stations, a mean-field asymptotic behavior is obtained for the dynamics of the 
distribution of vehicles over the different sites. Using the same method, [9] extends 
the study to systems where inhomogeneity is modeled by clusters. In addition, 
different alternatives are investigated, as for example, avoiding empty (respectively 
full) stations when taking (respectively returning) a bike, or returning bikes at the 
station having the largest available space, between two randomly chosen stations. 

The present paper analyzes the impact of finite capacity on large inhomogeneous 
networks. The models considered thus involve real-world blocking and rerouting 
mechanisms, that we implement on the infinite capacity model of [6] and [10j . 
The crucial product-forms of the stationary states are proved for a class of state- 
dependent routings that extends the classical setting of (5j. Asymptotics are ob¬ 
tained as both the number of vehicles and the size of the network increase. Namely, 
it is proved that the finite dimensional distributions of the stations occupation pro¬ 
cess at stationarity converge to products of truncated geometric distributions. If 
the total number of bikes does not exceed some order fo magnitude, stations and 
routes are moreover asymptotically independent, the latter being approximately 
distributed as Poisson variables. The approach, at stationarity, is similar to [6] and 
[TOl . As compared to the dynamical point of view of [8] and [9], our results com¬ 
plete the picture of the equilibrium state there obtained for locally homogeneous 
systems. 
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Our asymptotics fit the frame of the so-called principle of equivalence of canon¬ 
ical and grand canonical ensembles, in the Statistical Mechanics terminology. The 
networks of interest are here in some “partially subcritical” phase, in the sense that 
part of the network - namely, the parking stations - is in normal traffic (or no¬ 
condensation) regime. This makes that a Local Limit Theorem can be used. This 
way of deriving the grand canonical approximation is standard and can be traced 
back to Khinchin m- Dobrushin and Tirozzi [4] have used it for analyzing a family 
of Gibbs measures on Z 2 . It has, since then, proved useful in many different con¬ 
texts, notably for exclusion and zero-range processes under thermodynamic limit 
m Appendix 2). Note that zero-range processes are identical to (homogeneous) 
closed Jackson networks, but due to homogeneity, condensation can only occur in 
specific attractive cases, which are not considered in [14] . Such supercritical phases 
are described in 12] and 11. The generalized Jackson networks with blocking con¬ 
sidered in Section [I] contain as particular cases the generalized exclusion processes 
analyzed in [14j . 

Contrary to the classical setting of the equivalence of ensembles - or its alter¬ 
native formulation as Gibbs conditioning principle ([3]) - the equilibrium states of 
general closed Jackson networks are non-symmetrical with respect to the nodes. 
Malyshev and Yakovlev ([IS]) have used analytical methods to address the prob¬ 
lem for classical (one-server nodes) closed Jackson networks, under existence of a 
limiting density of customers - as in [14] - and of a limiting profile, capturing inho¬ 
mogeneity, for the traffic parameters. The critical phenomenon of condensation is 
highlighted and both phases are studied. These results are extended in .6], using 
Local Limit Theorems, which is more flexible and does not require specification of 
a limiting density nor parameter profile. Note yet that this method does not lead 
to a complete and unified treatment of supercritical regimes, for which restrictive 
assumptions are needed. 

We use the same approach and for simplicity, restrict the analysis to the relevant 
dynamics in respect to the targeted vehicle-sharing applications. Though our inter¬ 
est focuses on finite capacity systems, the infinite capacity case - that could equally 
be derived from [6] - is included, as a reference and for completeness. Similarly, 
several standard results (as regards the Local Limit Theorem, or generalized Jack- 
son networks) are here stated, under specific assumptions adapted to the present 
context, for the sake of self-containment. 

The paper is organized as follows. Section [2] displays the Central Limit and 
Local Limit Theorems - for independent, non identically distributed, random vari¬ 
ables - that will be used in the sequel. Section [3] introduces the generalized Jackson 
processes of interest and states the product-forms of the associated stationary dis¬ 
tributions. In Section [4] the equivalence of ensembles is proved for two classes of 
networks (with respectively infinite and finite capacities). Applications to evaluat¬ 
ing the performance of bike-sharing systems is finally developed in Section [5] 

2. A Local Limit Theorem 

It is classical that the Central Limit Theorem holds for any sequence (Sn)n>i 
of sums of independent square integrable random variables 

J(N) 
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if the family (Xj } N)i<j<j(N) satisfies the following so-called Lyapunov condition: 

1 J(N) 

(2.1) 35 > 0 such that lim 2+g ^ E (| Xj t N — mj j jv| 2+ ' 5 ) = 0. 

N ^°° b N j =1 

Here and in the sequel, for all TV > 1, J(TV) is some integer such that J(N) > 1, 
and we use the following notations: 


J(N ) 

mj^N = E(TV^at) (1 < j < J(TV)), On = E(Sjv) = m j,N and 

3 =i 


J(N) 

&j, N = E([X,- iJV - rnp N ] 2 ) , b% = E ([S N - a N ] 2 ) = ^ a j,N- 

3=1 

This result is usually derived from the Lindeberg Central Limit Theorem, which 
states that under the following Lindeberg condition m , which is weaker than 
(12.11) . the characteristic function (or Fourier transform) of 6^ 1 (5 jv — a at) converges 
to that of the standard normal distribution (see [2|, or |7j): 


(2.2) Ve > 0, 


lim ^ VE ((Xj n - TOyjv) 2 

AT->00 bi T v J J ’ 

iv j — 1 




It is easily seen that, under the Lyapunov condition, the convergence is uniform 
on some sequence of intervals [— which length grows to infinity. This 

refinement of the Central Limit Theorem is stated below as Lyapunov Central Limit 
Theorem. It will be used for deriving the next Local Limit Theorem that is the key 
to the equivalence of ensembles. 

The proofs of both theorems, as well as those of the two next propositions are 
deferred to the Appendix. 


Let (X/ i jv)iv>i,i<j<j(iv) be a family of square integrable random variables such 
that, for all TV > 1, ..., Aj(at),at are independent. Then define Sn,cln and 

b at as above. 


Theorem 2.1. (Lyapunov Central Limit Theorem) 

If the Lyapunov condition en) is satisfied, then the following convergence holds, 


(2.3) 


lim sup 

|t|<Ajv 



— e 


-t 2 /2 


= 0 , 


for some sequence (An) of positive real numbers converging to infinity. 


All the models of interest in this paper will satisfy the Lyapunov condition with 
5 = 1. 


Deriving of the Local Limit Theorem from the Central Limit Theorem for integer 
valued random variables has become standard since the paper by Gnedenko Eli!, 
addressing the case of i.i.d. variables. We here give a version adapted to our context. 

Theorem 2.2. (Local Limit Theorem) 

Assume that the random variables X^n are integer X-valued and that the follow¬ 
ing conditions hold: 

(1) liniAr-Kx, bw = Too, 
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(2) there exists some sequence (An) of positive real numbers converging to in¬ 
finity such that (USD holds, 

(3) there exists some <fi £ L 1 (K.) such that for all N > 1 and all t £ [—7r,7r], 

|E(e itS ")| < cj>(b N t). 


Then the following is true 
(2-4) 


lim 

N—too 


sup 

fcez 


b N V2nP(S N = k) - exp 1 j 


= 0. 


We will now consider two examples of interest for our applications of Section [4] 
Conditions (i) of the following Propositions 12.11 and 12.21 will there mean that the 
networks are subcritical as concerns their one-server nodes. 


Example 1. 

Suppose here that, for all N, the set of indices {1,..., J(7V)} is the union of two 
disjoint sets Jfj and J'f, such that 

• for j £ the random variable Xjn is geometric with parameter pj ,n £ [0,1[, 

i.e. P (X jtN = k) = (1 - Pj,N)pj, N i for k e N, 

• for j £ J'fj, the random variable Xj t N is Poisson with parameter A j ,N > 0 , 

A fe 

i.e. P(X 7 - N = k) = e~ Xj ' N f ,JV , for k £ N. 

k\ 

The means and variances of the X^jv’s and S/v’s are then given by 

Pi.N r ■ p qrN ( Pj,N r ■ p qrN 

1 - Pi jv 101 3 -2 J ( 1 - 0 , IOr J -'l 


mj ' N \ A 3 P n’ N for j £ J 2 n ’ 


?In=\ 

1 X J,N 


for j £ Jf 


N > 


on = 


E Pj,N 

1 — Pi N 

jeJ N n, 


\' 


\' 


Pj,N 


X hN and b N = ^ n _ n . „\2 


\ ' 




j&JT ^ Pj ’ N jejf 


A, 


*j,N- 


Proposition 2.1. If the two following conditions are satisfied: 

(i) there exists some p < 1 such that for all N >1 and j £ pj t N < p, 

(ii) linuv^oo b N = +oo, 
then m holds. 

Remark 2.1. (*) and (ii) of Proposition HOI are not necessary conditions: The 
Local Limit Theorem can still hold beyond the subcritical phase. It can be proved 
for example that if the number of indices j £ such that pj,N = ma x^jn pj^N 
goes to infinity and lim inf max je j-w pyiv > 0, then equations (12.31) and (12.41) hold. 

Example 2. 

Suppose now that, for all N, the set {1,..., J(N)} is the union of two disjoint 
sets and such that 

• for j £ the random variable X^n has truncated geometric distribution 
with parameters pj,N > 0 and Cj t N £ N, i.e. for k £ {0,..., Cyjv}, 
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• for j £ the random variable X^n is Poisson with parameter Ayjv > 0 
as in example 1. 

The means and variances of Sn’s are here given by 

aj\r = 'y ] 171 j,N + y ' Ayjv and b 2 N = y ] a j,N + ^ ' Ayjv, 

j&J? jeJ 2 N jejf j&J 2 N 

- (Cj.JV + l) Pl 'i+c j N if Pj,N 7 ^ 1 
1 -Pj,N 

if Pj,N = 1 


where for j £ J™, 


Pi, N 


7 j,N = 


- I i- 


Pj,N 


Cj.iV 

2 


and 


J2 

a 0,N 


Pj,N 


(1 -Pi.NV (°3’ N 

C .7,-W(c.7,jV+2) 

12 


i) 2 


Pj,N 3 ' 


1 + c j,N \ 2 

l 1 P i,N ) 


if Pj,N ^ 1 
if Pj, N = 1 


Proposition 2.2. If the two following conditions are satisfied: 

(i) there exists some C < +oo such that for all N >1 and j £ Cyjv < C, 

(ii) lim^v-s-oo b N = +oo, 
then 422} holds. 

The proof of this proposition (see Appendix) requires uniform domination of the 
characteristic functions of truncated geometric variables. This is stated as Lemma 
12.21 which itself relies on the following lemma, due to Gnedenko El p. 192-3]. 

Lemma 2.1. (Gnedenko) Let X he an integer-valued random variable. Denote 
p n = P(X = n ) for n £ Z and 

_ P 21 P 21+1 

lez P21 + P2i+i ’ 

with the convention that p ^+p'] hl 1 =0 if P 21 = P 21+1 = 0. Then, the following 
inequality holds 

(2.5) | R(e itx ) | < e~^ st2 for \t\<n. 


Now for any p > 0 and c £ N, let m PjC and a 2 c denote the mean and variance of 
a truncated geometric random variable X PtC with parameters p and c, and let s PiC 
denote the associated parameter as defined in Lemma 12. II The following lemma is 
easily derived from Lemma 12.11 


Lemma 2.2. For any C > 1, there exists some n > 0 such that for any p > 0 and 
any c £ {0,..., C}, 

| E(e itXp ’ c ) | < e -^ K <<= t2 f or \t\<TT. 

Proof. It is enough to consider the case when 0 < p < 1, since a 2 p _ x c = a 2 c and 

| E(e itXp ’ c ) | = | E(e it(c- ' Y p -1 ’ c ) | = | E(e _i ^- 1 .<=) | 

due to the following elementary remark. 


Remark 2.2. For all c £ N and p £]0,+oo[, the random variable X PtC has the 
same distribution as c — X p -i c . As a consequence, the following symmetries hold 
for any c £ N and p > 0: 


TYlp— i ? c — ^ TYlp,c 



= a 


2 

p,c* 


and 
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This induces a duality property between occupied and empty rooms in the finite 
capacity systems to come, known as the particles-holes duality (see [14]). 


So, it is now enough to prove existence of some positive k such that 

s p,c > K (: r p,c 

for all c € {1,..., C} and 0 < p < 1. (The case c = 0 needs not be considered, 
since sq, p = <Jq p = 0 for all p.) 

It is clear that for fixed c > 1, both mappings p K > s PjC and p a pc are positive 
valued and continuous on the interval ]0,1]. It is then enough to prove that, for 
c = 1,..., C, 

liminf s PiC /a 2 0 c > 0. 

p—>0 ' p ’ c 

Now from the definition of s P)C , one gets (using p < 1 for the second inequality) 


*c= E 


n 4i+l 


> 


\l =0 
P 


, 0<22<c-l 


p2l p2l + l 


E p 21 E^ 


> 


, 0<2Z<c-l 


G=0 


E^ 

U=0 / 


Together with a p c = — (c+ l) 2 (1 ^ pC+ i^ if p < 1, we obtain 


lim inf -Pff- > U m inf 
p->o o z pc p^> o 




G=0 


(1 -P? 




(i - p^y 


1 

2 


□ 


3. Generalized closed Jackson networks 

In accordance with the terminology of Serfozo nn, a standard Jackson network 
consists of a finite number of nodes at which customers are successively served, 
according to possibly different procedures. Nodes operate independently, which 
means that the rate at which customers leave a node only depends on the current 
occupation of this node. When a customer completes his service at some node, 
he instantly moves to another node for another service. New nodes are chosen at 
random according to a fixed routing matrix. All service processes and routings are 
independent. For a closed Jackson network, which is the case considered in this 
paper, customers stay forever in the system and there are no external arrivals. 

In this description all nodes have unlimited capacity. This dynamics can be 
generalized to a system including nodes with finite capacities. But the transition 
rule needs then be state dependent (constrained by the current location of the free 
space). It can no longer be given by some fixed matrix. Different blocking or 
rerouting policies can be considered. 

3.1. Infinite capacity model. First consider the case when all node capacities 
are infinite. Denote by N the number of nodes, by M the fixed number of customers 
and by P = (p*j)i<*,j<jv the routing matrix (that is, p,j is the probability that a 
customer leaving node i moves to node j). P is assumed irreducible, which ensures 
uniqueness of its associated invariant distribution 6 = (&i,, 8 n) on {1,..., N}. 
For each node i, the departure rate of customers from node i is gi(k) when k 
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customers are present at i. Here gi is some function defined on N such that gi{ 0) = 0 
and gi{k) > 0 for k > 1. 

The node-occupation process is Markov with state space 

rij =m\ , 

and transitions 


f N 

Sn,m = < (ni,..., un) £ such that ^ 

l i=i 


n —> n — ei + e 3 at rate q(n , n — + e 3 ) = gi(rii)pij , 


for n £ Sn,m and 1 < i,j < N. Here d denotes the i th unit vector. 

A well known and remarkable feature of this class of processes is that the joint 
stationary distribution of the queues lengths £i,..., £jv is explicitly known and has 
product form, given by 


^ N Q n o 

(3.1) P(£i = m,... ,£n = un ) = tz TT 7—ITT—T for eSjv.Ar, 

Z fJi w')K) 


where Z is a normalizing constant and for j = 1,..., N, the function (gj!) is defined 
on N by 

n 

fe ! )(0) = 1 and (ffj!)(n) = g 3 (k) for n> 1. 

fc=i 


3.2. Finite capacity models. If in the contrary some nodes have finite capacities, 
denote ct the capacity of node i (0 < c* < +oo for i = 1,..., N) and c the vector 
(ci,..., c/v). The state space is now changed for 

f N | 

Sn,m = ^ ( n ij • ■ •) n iv) 6 rij = M and rij < c 3 for c 3 < +oo > , 

l i =1 J 

N 

which is non empty provided that Cj > M . Several dynamics are considered. 

1=i 

1 - Model with blocking. 

The dynamics is the same as in the previous infinite capacity case, except that 
when node j is at capacity, any move from some node i to node j is cancelled (that 
is, the customer involved renews service at node i). 

In other terms, the transitions are given by 

n — >n-ei + e 3 at rate q(n,n - e* + e,-) = gi(ni)pijl nj<CJ , 
for n € S“ff M and 1 < i,j < N. 

The infinite capacity model can be seen as a particular case of this finite capacity 
model with blocking. But the explicit formula (13.11) for the invariant distribution 
does not extend, in general, to the finite capacity model with blocking. It is the 
case, yet, if the routing matrix P is reversible with respect to its invariant measure 
0, that is, if 

(3.2) OiPij = djpji for 1 < i, j < N. 

If (13.21) is satisfied, then (13.11) still holds here for (m,..., njv) £ S% M (see [16]). 


2 - Model with blocking and rerouting. 
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This is a particular case of a dynamics introduced by Economous and Fakinos 
[5]. Here, any customer routed - according to the matrix P - to some saturated 
node will move instantly (i.e., at infinite speed) across the system, still according 
to P, up to find some non saturated node where to settle for a new service. Note 
that this customer will certainly find such a node, possibly the actual node he has 
just leaved (where one unit of capacity is thus available). Of course, in this case, 
no transition occurs. 

The corresponding process has the following transitions and rates 

(3.3) n —> n — ei + ej at rate q(n, n — ej + ej) = gi(rii)p*j(n), 

for n £ M and 1 < i ^ j < N, where p*j ( n ) denotes the probability that a 
Markov chain with transition matrix P initiated at i enters the set Aj (n) at j , for 

(Ig. f 

Ai{n ) = {j, 1 < j < N, j =i or n 7 < c,}. 

Contrary to the simple blocking case where P has to be assumed reversible, here 
the product form ED holds without any restriction on the routing matrix P. This 
dynamics thus appears as the appropriate generalization of the standard - infinite 
capacity - Jackson dynamics (in which case, p*j(n) = ptj for all n,i,j). 

Theorem 3.1. (Economou and Fakinos). For any irreducible routing matrix P, 
the Markov process with state space , M and transition rates given by (13.311 has 
invariant distribution given by (ED for (m ,..., njsr ) £ Sfj M . 


3 - Other state-dependent routings. 

The scheme introduced in [5] is slightly more general than the blocking and 
rerouting dynamics just described. In the general setting of [5j, the set Ai(n) - of 
nodes that are below capacity when some customer has just left node i - is splitted 
into disjoints blocks. A partition is thus associated with each configuration n — ej 
of M — 1 customers over the N nodes. The just served customer at i still explores 
the network at infinite speed according to some Markov chain with transition P, 
but settles at the first visited node that belongs to the same block B n - ei ( i ) as i. 
The same product form stationary distribution then holds, given by (ED- 

This result can be generalized to a larger class of state-dependent routings. We 
consider the following transitions and rates: 

(3.4) n — >n-ei + ej at rate q(n, n- e, + ej) = gi{ni)p i3 (n — e») l ni>0 , nj<Cj , 

for n € Sfi M and 1 < i j < N, where a transition matrix P(m) = (pij 
on the set 

A(m) d = {j, l<j<N, mj < cj} 
is associated with each m £ 5^ M -i- 

Theorem 3.2. Assume that some positive vector 9 = (9 1 ,... ,9n) is such that 
Vm £ Vj £ {1,..., N} 0iPi j {m)=0 j . 

i£A(m) 

then the Markov process with state space M and transition rates given by (ED 
has invariant distribution given by (EJ for (m,..., un) £ Sfj M . 
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The dynamics considered in [5j are of the form described in (13.41) , with line i of 
matrix P(m) then given, for all m £ S^ M _ x and i € A{m), by the distribution of 
the return point to B n _ ei ( i ) of a Markov chain with transition matrix P initiated 
at i. The condition of Theorem l3.2l is then satisfied for the invariant vector 9 of P, 
due to the well-known following result: For any ergodic Markov chain on {1,, IV} 
with invariant vector (0i)i<i<jv, and for any subset B of {1,..., N }, the embedded 
Markov chain at times of visits to B has invariant vector (Oi)i£B- 


Proof. We check that the global balance equations 

(3.5) ^2 gi(ni)pij(n- ef) = ^ gj(nj + 1 )pji(n - a) ^ f\ + &3 \ 

rii>0,nj<.Cj n,i>0,nj<Cj 

are satisfied for n £ where n(n) = Z ” 1 IljLi /(9j■)( n j)- 

Recall that gifrii) = 0 if rii = 0, so that the left-hand side of (13.51) is equal to 


N 


^ 9 i{ n i ) ( X ~Pii( n ~ e i)) 


while, using 7r(n — e* + ej)/w(n) = Ojgifrii) / (0igj(rij + 1 )), the right-hand side 
rewrites 

5Z 7T9i{ni)Pji(n - e i)=Yl 55 e iPji( n - e 0 


i<4^j<iv 

n 3<Cj 


6i 

i—l rij <Cj 


Now equation (13.51) results from the following equality: for all n £ Sfj M and 
i £ {1,..., N} such that m > 0, 

Oi (1 - pain. - et)) = ^ 9jPji(n-ei), 


or equivalently (adding 9ipn(n — ei) to both members) 

9i — ^ ( ®jPji[ri Ci), 

iGA(n-ei) 

which is satisfied by assumption on vector 9. 

□ 


Remark 3.1. It will next be crucial to note that due to the particular form of 
the constraint Ylf=i n j = the expression in (13.11) . for (ni,... ,hn) in Sn,m or 
Sff M , is unchanged if the invariant vector 9 is replaced by some proportional vector 
7 9 = (g9 1 ,... ,j9n), where 7 > 0 is arbitrary. 


4. Equivalence of ensembles for generalized Jackson networks 


The product formula 63 can be interpreted as follows: the invariant joint 
distribution of the N queues is the distribution of N independent random variables 
conditioned to the constraint that their sum is M. 

Moreover, for arbitrary positive 7 , remark [XT] allows one to choose the individual 
distributions of these independent variables, denoted gj ,..., gjj, as given by 


(4.1) 


P(?yJ = n) 


1 ) n 

z Ai) (ffjON 


for 0 < n < Cj + 1, j = 1,..., N, 
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provided that 7 belongs to the domain of convergence of each power series 


(4.2) 


ZAl) = E 


91 


(9j !)(«) 


for j = l,...,N. 


n—0 


For each such 7 , JO rewrites: for (m,... ,un) G Sn,m, 
(4.3) P(£i =ni,...,6v=rajv) = 



N \ 

II 

..,r)]f = n N \ E^ 7 = M I 


In the physical terminology, the distribution in (13.111 is known as the canonical 
ensemble , while the product of the N distributions in mu is the grand canonical 
ensemble. The factor 7 is the so-called chemical potential. 

If it holds, the principle of equivalence of canonical and grand canonical ensem¬ 
bles tells that if 7 is rightly chosen, namely, such that for rfj’’ s as in (HU), 


(4.4) 


E 


N 

E 5 

ki =1 


= M 


then the finite dimensional marginals of the distribution m are well approxi¬ 
mated, for large N and M , by products of distributions in (HID. In other terms, 
for N and M large, one can forget the conditioning r/f + ■ ■ ■ + rjJ^ = M provided 
that 7 is rightly tuned, so that the total mean of the N free variables rfl’’ s is M. 

The following lemma proves strict monotonicity of the left-hand side in (14.41) . 
ensuring uniqueness of 7 solving this equation. Existence of such a 7 holds, by 
continuity, provided that for at least one j (such that, necessarily, Zj( 7 ) has radius 
of convergence 7 *) 

lim^ E ( 77 J) = +00 

where 7 * is the minimal radius of convergence of all the power series Zj( 7 ). 

Lemma 4.1. Let (X 7 )o <7<7 » be a family of N-valued random variables with dis¬ 
tributions 

P(X 7 = n) = 7 "f^ for n G N, 

Z (V 

where (f> is any nonnegative function on N that is non-zero for at least two points, 

OO 

z( a) = E^ n )’ 

n —0 

and 7 * is the radius of convergence of the series Z(j). Then E(X 7 ) < +00 for all 
7 G]0, 7 *[ and the mapping 7 G]0, 7 *[ 1 —> E(X 7 ) is increasing. 

Proof. Both series ^i=o n 'y n< t , ( n ) and E^Eo n 4>(n) have same radius of conver¬ 
gence 7 * as Z(j), and E(X 7 ) is given for all 7 G [0, 7 *[ by 

E(Eo 


E(X 7 ) = 


z ( 7) 


Differentiating with respect to 7 G ]0, 7 *[ gives 


<9E(X 7 ) 

~ 7 


E“o"V0(n) 

Z{ 7) 


E”o 

Z( 7 ) 


= 7 . VarX 7 > 0 
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from assumption that <f>{n) > 0 for at least two values of n, ensuring that all 
variables X 7 for 0 < 7 < 7 * are non a. s. constant. 

□ 

From now on, we will denote 771 ,..., 77 jv the variables r/J,..., rp N associated with 
the 7 solving COD , if it exists. They will be referred as the “free variables”. 


The Local limit Theorem is a classical tool for proving that equivalence of en¬ 
sembles holds m, a, m, M)- One can assume without loss of generality that 
the finite dimensional distribution of interest is that of (£ 1 ,..., £k) for some K > 1, 
and write for any ( 711 ,..., uk) G N a ": 



N \ 

II 

3 

■ jVk = n K | 55 7?J = m \ 

3 = 1 / 


= P (?7 i = 77.1, ... ,rj K = n K ) ■ 


p(e 7 =i^ = m) 


where the last equality results from independence of the 77 ^ ’s. 

If the Local Limit Theorem holds for both families of variables ( r]j , 1 < j < N) 
and (r]j , K + 1 < j < N), as N goes to infinity, then since by choice of 7 , 

' N 

E(? 7 i) + • • • + E(? 7 jv) = M, one gets, informally: 
and 



1 


N 


K 


exp 


% = m -5e 




&Ar\/27r 

(£f=i(Efai)-»i)) 


0=^+1 


j =1 




where b 2 N = cr|, denoting er ? - the standard deviation of rjj. 

,k 2 


Now if both Ej=i an d (EjLiQEfe) — n j )) 2 are negligible with respect to b 2 N , 
the right-hand side of the last approximation is close to (biyy/ZK)- 1 and one gets 
the expected equivalence 

P (£1 = ni, ■ ■ ■, £k = n K ) ~ P (771 = nr,..., r) N = n N ). 


To make this formal, we must consider a sequence of networks. In what follows, 
this sequence will be indexed by the number N of nodes, thus adding a subscript 
N to all parameters and variables. 


Remark that if the Local Limit Theorem applies to the complete vector of free 
variables (771,..., 77^), it provides an equivalent for the partition function Z in (13.11) . 
as function of (the implicit) 7 . Indeed, using equations (13.11) . (14.11) and (14.31) . 

( N \ ( N \ N 

n 7 ) p (=m »( 7 M b N v^r i n Zih). 
i= 1 J \i=l / i= 1 


For application to vehicle-sharing systems which is the object of Section [5l we 
need consider only two types of nodes: single server nodes (with finite or infinite 
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capacity) and infinite server nodes with infinite capacity. For a single server node 
j. when n customers are present, the departure rate from j is gj(n) given by 

gj(n) = if n ^ 0 and gj( 0) = 0, 

while for an infinite server node, gj is given by 

gj(n) = /JjTi for all n£N, 

where in both cases, /ij > 0 is the parameter of the exponential services at j. 
Equation ED then tells that the variable rjj is 

• geometric with parameter 7 t L J 1 0j, truncated to {0 ,..., Cj} if Cj < 00, for 
a single server node j, 

• Poisson with parameter 7 pj 1 0j if node j has infinitely many servers. 

We will first consider a model where all nodes have infinite capacity, and then a 
model where single server nodes have finite (uniformly bounded) capacity. The 
equivalence of ensembles for those two models will respectively be derived from the 
local limit theorems stated in Propositions 12.11 and 12.21 of Section [2] 


4.1. Networks with infinite capacity nodes. We here consider a sequence of 
standard closed Jackson networks: all nodes have infinite capacities. The network 
numbered N has N nodes (labelled by 1,...,7V) and Mjv customers. Nodes are 
divided into two disjoint sets and ■ Those in are single server nodes 
and those in J.± are infinite servers nodes. To avoid trivialities, the set J{* is 
assumed to be non empty. 

For each network, we add a subscript N to all quantities already defined for a 
general Jackson network, thus denoting 

• Pn = {pij,N)i<i,j<N the (irreducible) routing matrix, 

• On = (6 l i j jv, ..., &n,n) the invariant probability vector associated to Pn, 

• g>j,N the rate of exponential services at node j for 1 < j < N, 

• (& a) 1 <j<N the stationary node-occupation random vector. 

The product formula (13. If) here writes 

(4.6) P(£i,iv = nii • • •, Cv ,n = tin) = -5— 

ZjN 

where (m,..., njv) € Sn,m n and 

rj,N = ^ N 0 hN for 

The parameter rj n is usually referred to as the utilization of node j. From 
Remark 1331 one knows that those utilizations need only be defined up to a con¬ 
stant. For this reason, some authors normalize the r^jv’s in such a way that for 
each N, the maximum max{771 < j < N} is equal to 1. This normalization 
is crucial in [15], where it is fully part of the assumption that the empirical dis¬ 
tribution iW 1 £ i= i 6 rj N converges to some probability measure. Here, it will not 
be necessary, since less stringent conditions are needed: The factor 7 will in some 
sense absorb the normalization factor. 


n 

jeJ? 


j,N 


n 

2 n 


'j,N 


1 <j<N. 
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For all N, let qjv £ [0,1/maxjr^Ar, j £ }[ be the (unique, from lemma l4~Tl) 

solution to equation (14.41) . here given by 


(4.7) 


E 

ieJ? 


iNTj'N 
1 ~ lNTj,N 


E 7 Nr 3,N 
j£J 2 N 


m n . 


(The left-hand side of this equation clearly goes from 0 to +00 as 7^1 varies from 0 
to 1/maxjr^jv, j £ ensuring existence of the solution.) 

The free variables ( i)j,N)i<j<N associated to (£j,iv)i<j<jv are independent and 

• for j £ rjj t N is geometric with parameter 7 Nfj,N, 

• for j £ J 2 1 Vj,N is Poisson with parameter ''/NT'j.N, 

so that Ylf=i Vj,N has mean Mn (from equations (14.71) 1 and variance 


(4.8) 


b% 


E 

jeJ? 


lNrj,N 
(1 - 7lvTj,Af) 2 


E 'VK r 3,N- 


Theorem 4.1. Assume that the following conditions are satisfied 

( 1 ) limsup I 7 jv max{r 3 ,at} < 1 , 

iv —^00 V jejf J 

(2) lim 6 at = + 00 , 

JV-> 00 

(3) K > 1 is such that lim bfj 1 V jN'fj n = 0 . 

N—>oo ^— J 

Then for any (ni ,..., uk ) £ , 

f 49) Jim p (6 ,N = n 1 ,...,Z K ,N = n K ) 

1 j N—>oo P( m , N = m)... P( Vk , n = n K ) 

Proof. We use (|4.5I) (here with a second subscript N for all random variables). 
From equations (14.71) and (14.81) . Xq=i Vi has mean M N and variance b 2 N . 

Condition (1) together with inequalities 7 jv < 1/ma x^jn rj t N for all N imply 
existence of some p < 1 such that 


VW, Vj £ J”, lN r h N<P- 


Condition (2) then allows to apply Proposition 12.II to both families (' r)j,N)i<j<N 
and ( r]K+j,N)i<j<N-K (with respectively, J(N) = N and J(N) = N — K). This 
is immediate for the first family. As for the second family, it results from 


Var 


( E 

\j=K +1 J 


— b% — 


TN r j,N 
(1 — 'YNUi n) 2 

jej? n{i,. v 3 ’ ' 


E 7 Nr jt N, 


that Var \ -K+i Vj,Nj ~ b 2 N as N goes to infinity, since both sums in the right- 

hand side are negligible with respect to 6 ^: indeed the first sum is bounded above 
by Kp/{1 — p) 2 , and the second one is negligible with respect to 6 /v, hence to b 2 N , 
from assumption (3). 
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Then using (12.41) , 6 jvvP (j2^=i Vj, n — converges to 1 and 


\ 


K 

/ JV 

K \ 

2 ^-EEv) p 

y opN = m n 

E "v 

i=i 

\j=K+l 

1=1 J 


exp 


\ 


(r,f=i( m j,N - rij)) 

2 - Ef =1 °In) 


converges to 0 for any m,..., nx £ N, as N goes to infinity, where for 1 < j < K, 
m jt N = E( Vj,N ) and a 2 N = Var(r] jtN ). 

Now, on one hand, b% - a j,N = V ar (EyE+i %iv) ~ b 2 N as N goes to 
infinity, as just shown. On the other hand, the argument of the exponential in the 
last expression goes to zero since 

Y K K 

-^-J2\mj,N-n j \<bJ f 1 J2 n j+Kb] v 1 p/(l-p) + bJ i 1 ^ lNr jlN 
N i=i j= i jeJ^n{i,...,K} 

which goes to zero as N tends to infinity, using assumptions (2) and (3). Hence 

( N K \ 

lim b]srV2n P E = m n~Y 1 rij = 1 . 

\j=K+l 1=1 / 

Equation (14.91) follows, using (14.51) together with independence of the r]j : jv’s. □ 

Remark 4.1. (i) Note that, since M jv < b 2 N for all N, as can be seen from (14.71) 
and flUD , condition (2) of Theorem \f. 1\ is satisfied if limyv—>oo -Vv = oo. 

(ii) Assumption (3) is trivially satisfied for {1,... ,K} C . 


4.2. Networks with both finite and infinite capacity nodes. The networks 
now considered differ from the previous ones on two main points. For all N, 

(1) for j = 1 ,..., N, node j has capacity Cj t jv, which is now finite for j £ J\ , 
but remains infinite for j £ jJf 1 ; we will denote c/v = {cj,N■, 1 < j < N), 

(2) both Jj ,v and J™ are assumed to be non-empty. 

Any of the dynamics described in the previous section can be considered: block¬ 
ing, blocking-rerouting or general state-dependent routing satisfying the setting of 
theorem In the simple blocking case, we assume that for all N, the routing 
matrix is reversible with respect to its invariant distribution. Thus, in any of these 
situations, the queue-length process has stationary distribution given by the same 
formula (14.61) as in the infinite capacity case, here with state space Sff MN . 

We keep the same notations as before: ry yv are the (non-normalized) utilizations, 
Pj.N the exponential service parameters, On the invariant probability vector, either 
of the fixed routing matrix Pyv (for blocking or blocking/rerouting procedures), or 
common to all state-dependent routings, in the setting of theorem 13.21 

For all N, ~/n is now uniquely defined in [0, + 00 [ by 


E 


Snto n {lNTj, N ) 

££= 0(7 N^r 


E ^rj, N 

j£J 2 N 


M n , 


(4.10) 
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or equivalently by 



lNTj,N 
- lNTj,N 


( c j,N + 1 ) 


( 7 Nr JtN ) Cj ’ N+1 

1 - ( 7 ivrj,Ar ) C3 '' JV+1 


'YN'f'j,N = Mn 


with the following abuse: term j of the first sum has to be replaced by Cyjv /2 if 
lNTj,N = I- 

Existence and uniqueness of yjv again result from lemma 14.11 and from the fact 
that the second term in the left-hand side of equation (14.101) increases to infinity 
with 7 n (recall that yf 0 ). 

The free variables (pj,N)i<j<N associated to (£,j t N)i<j<N are now such that 

• for j £ , r)j t N is geometric with parameter 7 Nfj,N truncated to [ 0 , Cyjv], 

• for j £ J 2 i Vj,N is Poisson with parameter 7 Ni"j,N- 

Denoting by and a 2 N the mean and variance of i)j t N (for 1 < j < iV), equation 

(14.101) rewrites Mn = 1 ^j,N, while b 2 N = a j n or equivalently 


b N ~ 


E 




(1 - 7 a Tr jtN y 


— (Cj,N + 1)^ 


(■ Int^nY 


i, JV + 1 


(1 - ( 7 N rj, N ) Cj ’ N+1 y 


+ 


E ^ Tr i’N- 


Here again, term j of the first sum must be replaced by Cj t N{cj,N + 2)/12 when 
lNTj,N = I- 


Theorem 4.2. If the following conditions are satisfied 

(1) there exists some C < +00 such that for all N > 1 and j £ Cj t N < C , 

( 2 ) limjv-^oc b N = +00, 

(3) K > 1 is such that lim bf^ 1 Int^n = 0 . 

then (POl) holds for all (rii ,..., Uk) £ N A ". 

Proof. The proof is the same as for Theorem l4.ll here using Proposition 12.21 Con¬ 
dition (1) ensures that all to^jv and ct^jv for j £ J-fi are bounded by some constant 
(not depending on N), namely, by C 2 (since the corresponding rjj jv’s have values 
in {0,..., C}). ’ □ 


Note that ( ii ) of Remark ro still holds. But this is not the case for ( i ): if 
m PiC and a 2 c denote the mean and variance of some truncated geometric random 
variable X PtC with parameters p and c, the ratio a PtC /m PtC approaches zero as p 
goes to zero or to infinity, for fixed c > 1 . 


5. Application to bike-sharing systems 

The purpose of this section is to derive practical results on the performance of 
bike-sharing systems such as the Velib’ network in Paris. Here, a (large) set of bikes 
are distributed over a (large) number of docking stations and offered for use to the 
population of the city. Any user can take a bike at some station, and then ride to 
another station where he returns the bike. The payment process requires that each 
bike be locked to a terminal, so that each station can accommodate only a given 
number of bikes. 
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5.1. Infinite capacity approach. It is easily seen that such a system can be 
modeled - at least as long as all stations remain below capacity - by a closed 
Jackson network with a fixed number of customers (here the bikes), and both single 
server nodes (the stations) and infinite servers nodes (the routes from one station 
to another). Indeed, assuming that the arrival process of users at any station a 
is Poisson with parameter p a , then bikes leave station a at rate /i a (here, inter¬ 
arrival intervals between users act as service durations for the waiting bikes). Next, 
assuming that ride durations are independent and exponentially distributed with 
parameter p r for route r, the number of bikes on route r decreases from n to n — 1 
at rate p r n. 

These remarks lead to a model with infinite node capacities proposed in [6j and 
[ID] , that is considered now. Denote by Jx the set of stations, with J\ = \Jx\-. and 
J 2 C {[ij] : i,j £ Jx), with J 2 = \ Ji\ < Ji(Ji — 1), the set of possible routes: 
Here [ij] denotes the route from i to j. It is assumed that the bike moves obey to 
some statistics which is constant in time: that is, any user taking a bike at station i 
has probability qxj to put it back at station j. Stations here have unlimited parking 
capacity, as well as routes can accommodate as many riders as necessary. Then 
bikes/customers alternatively occupy nodes in Jx and in J 2 and move according to 
the following routing matrix P on J\ U J 2 : 

(5.1) Vi,j£Jx, Pi[ij] = and = 1 ; p ab = 0 for other a, b £ Jx U J 2 . 

Here J 2 is assumed to include all [ij] such that q^ > 0, so that (15.11) well defines a 
routing process on J\£)J 2 . All services (inter-arrival times between users at stations 
or trip durations) are independent, exponentially distributed with parameter p a at 
node a (a = j or [ij] for some i,j £ Jx) and independent of the routing processes. 
This model thus constitutes a standard closed Jackson network with M customers 
and N = Jy + J 2 nodes. 

Note that P is irreducible if and only if the following two conditions are satisfied: 

(1) the matrix Q = ( Qij)i,jeJi is irreducible, 

(2) Vi, j £ Jx, [ij] q[ij] > 0. 

We assume that this is the case. Then, denoting by v = (i'j)jeJi the unique 
Q- invariant distribution on J \, the routing matrix P has invariant distribution 
0 = (0 a )aeJiUj 2 given by 

(5.2) 0j = for j € Jx and %] = \for [ij] £ J 2 . 


The process so defined fits the frame of paragraph 13.11 of Section [3] and has 
stationary state given by 


(5.3) 


- "1.£.v - "x) - 7 n r 7 n 


jeJi 


[b']eJ 2 


1 M 

n [ij] ' 


where (nx, • • •, un) £ Sn,m and 

Vi 


for j £ Jx, 


2 p,j 


and for [ij] £ J 2 , = 


Vjqij 
’ ^[ij] 


This model does not wholly convey the complexity of real bike-sharing systems, 
since it ignores the blocking mechanism at saturated stations (which is the main 
problem for these networks). Still, it describes the system up to the first blocking 
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time. In [TO], asymptotic results are provided for M large, as N is fixed. Here 
Theorem 14.11 is used to describe this network as both M and N are large. 

The practical problem addressed is the following. For a network with given 
stations and routes, how many bikes and parking places should be offered ? It can 
be assumed that from observation-based estimations, the demand of the users is 
well known: that is, the frequencies at which users arrive at the different stations 
{fij' s), the popularities of the different routes (g^-’s) and their mean ride durations 
(p.jjjj’s). In other terms, utilizations r a = are known. Theorem 14.11 then 

helps evaluating, through approximation by geometric variables, the probability 
that some station is empty, or exceeds some level of occupation, indicating how to 
adapt the stations capacities and/to the number M of bikes offered. 

In order to derive results from Theorem 14.11 we again consider a sequence of 
networks indexed by N, and as in section [2 add a subscript N to all parameters 
defined above, rewriting equation (15.31) as equation (14.61) . 

A simple way to ensure that condition (1) of Theorem 14.11 is satisfied is to choose 
some small positive <5 and set 


(5.4) 7 N 


1-6 

maxj e jN r yjv 


and Mn 


E 

jeJ? 


lNrj,N 
1 - 7iV r j,AT 


E Tivr[y],jv, 

MeJ 2 N 


so that 7 jv solves equation (14.71) . Theorem 14.11 then has the following corollary. 

Corollary 5.1. Let 5 £ [0,1] be fixed and for all N, define jn and Mn by (15.41) . 
(i) If all rj t N for j £ Ji are of the same order, that is if 


mm jeJ if rj, N 

lim inf- J > 0 . 

iv->-oo maxj^jN rj t N 


then as N goes to infinity, the stationary queue-lengths at stations get asymptot¬ 
ically independent, and approximately geometric with parameters jN^j ,N = (1 - 
S) r j>N / max ieJ N r i>N (j £ J^). 

(ii) If moreover, 

max^jgj-jv r^jv 

lim — -- ---= 0 , 

N ^°° V Jl max t6 j» rk,N 

then at stationarity, both stations and routes become asymptotically independent, 
and the level of occupation of route [ij] is approximately Poisson with parameter 
lNT[ij] } N = (1 - <>) r [ij],N/ m ax feSi7i iV Tk,N- 

Remark 5.1. This corollary applies in particular under the following natural set 
of assumptions, for which 7 jv and Mn are both of the order of : 

• all rj t N for j £ Jare of the order of 1/ , 

• all r\ij] t N for j £ are of the order of 1 / i 

(as for example if all pj t N and P[ij],N ore of the order of 1, v^n of the order of 
1/Ji and v%,N<lij,N the order ofl/J % . Recall that JT v^n — "Yhhj] u i,NQ[ij],N = !)■ 
Indeed, y/j?/Jg < 1 /y/J? -+ 0 as N — > 00, since < J 2 by irreductibility of 

(q? j ),andN = J» + J?<(jn 2 . 
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Proof, (i) Due to (ii) of Remark 14.11 we need only check that condition (2) of 
Theorem 14.II is satisfied (condition (1) holds by choice of yjv). But here 


E 

jejf 


b 2 N > M n > V 7 Nfj^N = (1 — S) 


Jljej™ r i,N 
ma x jeJl r jtN 


>(1 ~S)Ji 




max jv 


which goes to infinity with N, since > y/N and lim inf w- 


min JV r jN 

— j£Jl > 0. 

™ jeJ iv r } , N 

(ii) We need to prove that max^jgjv T-W^ul.-v/fr/v S oes zer0 as ^ goes to 
infinity, in order that (3) of Theorem 14.11 holds. But as just shown, b 2 N > cJ± for 
some positive c and N sufficiently large. Hence, for large N, 

V[*j] G J2 1 lNr[ij], N /b N = (! - < 5 )- rhl] ' N ' 1 S 


b N m.ax keJ N r k ,N 
which maximum over [ij] goes to zero by assumption. 


< 


r m,N 


y/cJ f max kejf r k,N 


□ 


As a practical consequence, under condition (i) of Corollary 5.1, the probability 

• for station j to be empty is approximately 1 — (1 — S)rj t N / max j£ jx 17 jv, 
which is at least S , and exactly <5 at the most loaded station, 

• that station j has more than n vehicles is [(1 — '5) r yv/maxjgjj rj ) jv]” +1 . 
Capacity of station j can then be fixed as the smallest n such that this probability 
of exceeding level n is less than some given e. 

In the ideal situation when all utilizations ryjv for j € l7| ,v are equal, the above 
probabilities are the same for all stations. The probability for any station to be 
empty, 6 , is determined by the total number Mjy of bikes, given by equation (15.41) . 
And the smallest capacity that sets the blocking probability (at any station) below 
some given level e is given by the integer part of loge/log(l — 6). 


5.2. Finite capacity approach. We now directly consider networks with finite 
parking capacity (and still infinite capacity routes). Remark first that the routing 
matrix of the model of m, described in Section l5Tl is not reversible (p[ij]i = 0 for 
i ^ j while PiUj] = qij )■ Thus, the derived model with finite capacity stations and 
simple blocking procedure (as defined in Section liOl) is not tractable, its station¬ 
ary distribution being unknown. (Anyway, this procedure would not be realistic, 
making blocked users perform again the same ride as the one just completed). 

Still, two options are open for modeling such finite parking capacity networks. 
One is to consider the same routing as in [10] together with the blocking and rerout¬ 
ing policy described in Section [3] This choice turns out to be relevant for bike¬ 
sharing systems. Indeed, since transitions occur only from stations to routes and 
from routes to stations, any user blocked at his end-of-route station i will choose 
some route at random among the routes issuing from i, which well describes the re¬ 
ality. Note that in the standard blocking and rerouting procedure, the new route is 
chosen according to the same probabilities qij ’s as for new arrived users, which may 
seem unrealistic, or approximative (the q^s need then combine behaviors of new 
arrived and redirected users). We will discuss existence of other state-dependent 
routings that may be tractable (i.e. with product-form stationary state) and credi¬ 
ble for bike-sharing. 

The other option is to use pure blocking dynamics for a simpler model, in which 
all the routes are aggregated into one unique node, having infinite many servers 
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and infinite capacity. The transition probability from any station to this unique 
route is then equal to 1, while the transition from this route to any station j is 
given by some probability qj , called the popularity of station j. This basic model 
is a simple, easy to handle, approximation of the network in [10J . which does not 
make account of the detailed movements of the vehicles. It has been introduced 
and analyzed in [9]. Note that the associate routing matrix P is reversible, so that 
the product form (ED holds. Indeed, numbering the TV nodes (TV — 1 stations and 
one route) so that the set of stations is J\ = {1,..., TV — 1} and the unique route 
is TV, then P is given by 


(5.5) PjN = 1 and pNj = qj for j £ J\ ; pij = 0 otherwise, 
and the invariant probability 6 = (0 1; ..., 9 jy) by 

On = 7} anci 0 3 = ^ for j = 1,..., TV - 1, 

which clearly satisfies the reversibility condition (13.21) . 

Denote by Cj , for 1 < j < TV, the finite capacity of station j. The stationary 
state of the network is here given by 


P(£i = n\,... ,£jv = tin) 


1 

Zn 


r 


TIN 

N 


UN ! 


IV-1 

id; 


J=1 


for n £ such that n\ + • • • + njv = M and n 3 < Cj for 1 < j < N, where 
for j G {1,.. ., TV - 1}, rj = and r N = • 

Zi fjij A fJj 

Here, since there is one unique route, the simple blocking and the blocking and 
rerouting procedures are actually undistinguishable, both leading to the above prod¬ 
uct form. 


Remark 5.2. These simplified set of nodes and transition matrix can also be used 
in the infinite capacity case, instead of those of |10j. 


In modeling a bike-sharing system by any of these two models, it seems realistic 
to suppose that all the capacities c 3 stay bounded as N is large. Theorem 14.21 can 
then be applied. 

For the last model with one unique route, parameters /ii of services (inter-arrivals 
of users and durations of rides) should be assumed to be of the order of 1, while the 
qf s should be of the order of N~ 1 . One gets r 3 ’s of the order of TV -1 for 1 < j < N 
and rpf of the order of 1. 

Equation (14.101) here rewrites 


M = 



( c j + 1) 


(7 0) CJ+1 \ 

1 — (7 Tj) c i +1 ) 


7Hv 


(with the same abuse as in Section |TJ for terms j £ {1,... ,7V — 1} such that 
777 = 1). It can be used, as in the infinite capacity case, to fix M by choos¬ 
ing 7. The appropriate scale is here 7 = kN for some positive n, which gives 
asymptotically independent stations, with truncated-geometric approximate distri¬ 
butions with parameters n.Nrj and <7 (1 < j < N). Indeed, since b 2 N > 'yrw 
which is of the order of N, the approximation (14.91) of Theorem l4.2l is here valid for 
{!,..., iV} C Ji ={!,..., TV - 1}. 
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Note that nothing can be said about node N, since is of order y/N (jcn < 
b 2 N < 7r^v + (iV — 1) rnaxcj) so that bff'yr^ gets large with N. The unique route 
is thus not proved to be independent from the stations, nor to have some identified 
approximate distribution: In physical terms, this scaling of M (through 7 = kN) 
corresponds to a supercritical regime, with condensation at node N (that accomo¬ 
dates as many vehicles as the order of M). 

As regards stations: for 1 < j < N, 

• P(£j = 0) ~ 1/ En=o( K ^ r i) n > which decreases with respect to k 

• P(£j = Cj) « ( nNrj ) Ci / ^Zn=o( K ^ r j ) n i which increases with respect to n. 

If capacities are given, one can look for some k that gives a reasonable trade¬ 
off between the previous probabilities. These estimations can also help defining 
the capacities. Indeed for fixed k, the approximation 1/ Yjn=o( K ^ r j) n °f -P(£j = 
0) decreases as Cj increases. This is also the case for {uNr^Y 1 / E«=o( K -^ rr i)"’ 
that approximates P (£ 7 = c,), due to Remark 12.21 So Cj can be chosen as the 
smallest value that sets both quantities below some given level e. However the 
first quantity cannot be decreased below 1/E^o { K ^ r j) n (which is positive if 
uNrj < 1). Similarly, the second quantity cannot get below 1/E^Lo( K -^ r i) -n 
(positive if uNrj > 1). Here again, the ideal situation is when all 77 are equal. 
This indeed allows to select k such that nNr 3 = 1 for all j = 1,..., N — 1; so that 
the approximate values for P(£j = 0) and P(£j = Cj ) approach zero for large 
Cj’s. But real networks generally do not satisfy this condition, so that k can only 
be fixed such that nNrj = 1 for a group of stations with equal utilizations. Both 
probabilities - that a station is empty, or saturated - can then be arbitrarily reduced 
only for these given stations. No choice of k and Cj can be globally satisfactory. 

We will now consider models derived from that of [10], here limiting station 
capacities. We first investigate possible alternatives to the standard blocking and 
rerouting procedure - in which rerouting is ruled by the same matrix Q that rules 
choice of destination. The setup is the same as described in section l5Jl except that 
each node j in J\ has finite capacity Cj. The Jackson network dynamics is then 
modified, to avoid overflow of capacity, through the following transitions and rates: 
For n G <S^ M , where c = ( c a ) aeJlU j 2 , setting = oo for [ij] G J 2 , 

Ci -f- Rt rate Qij lni>o 

&[ij] “I - Cj rate ^[ij] lnj<cj 

C[ij] e [jk] at rate f^[ij] ij] C[%j]) lnj=cj ,n[ij]>o 

Here for all m G M _ 1 and j G Ji, a Markovian transition matrix W^\m) = 

( w ik ( TO )) on the set Ji\{j} is given. The following result is a consequence 

of Theorem 13.21 

Proposition 5.1. If for each j G J \, all matrices W^\m) for m G Sf, M _ 1 solve 
the following set of equations, with unknown variable W = (wik)i t keJi\{j} •' 

(5.7) Vk G J 1 such that k Y j, VjQjk = ^ v^Wik, 

ieJiUi} 

then, the process with state space S^ M and transitions defined in (15.61) has station¬ 
ary distribution given by (15.31) for (m,. .., niy) G SY M - 


(5.6) 


n 

n 


n 

n 
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Proof. It is easily seen that the rates in (15.611 have the form in (13.411 with 

for he N, gi(h) = /q l h>0 for i e Ji and g\ij]{h) = for [ij] G J 2 

and markovian transition matrix P(m) on the set A(in ) = {j e 3\,uij < Cj} U J 2 
given, for m G S% M _ V by 

Pj,[ij]( m ) = Qij for [*j] ^ Ji such that i G A{m ), 

P[ij],j( m ) = 1 for [ij] G J 2 such that j G A{m ), 

P[ij],[jk\(m) = u>ik(rn) for i,j,k G Ji such that [ij], [jk\ G J 2 and j ^ A(m), 

all other transitions having null rates. 

Moreover, it is straightforward to check that for all m G Sfj M _ 1; the restriction 
of vector 6 of (15.21) to the set A(m) is invariant with respect to P(m). This results 
from equations (15.71) together with invariance of v with respect to Q. 

The proposition then results from Theorem 13.21 

□ 

Equations <53 are clearly satisfied if = qjk for all i,j,k, which is the only 
solution such that depends only on ( j , k). In other terms, the above dynamics 
generalize the standard blocking and rerouting - for the Jackson network in jlOj 
- in the sense that redirection of users blocked at end of route [ij] may now take 
into account their original station i. Note that this excludes natural reroutings to 
stations in the neighborhood of j. Still, for each j , equations <53 have infinitely 
many solutions, as it appears by rewriting them as 

qjk = ^2 hi iW ik or else q :j . = qj, W, 

where Q = iffij) = (yiqij/i'j) is the time-reversed matrix of Q under its invariant 
vector v, and qj,, qj, denote the j-th lines of Q and Q. No general explicit solution is 
available, save for = (qjk)i^j,k^j- But in the special case when Q is reversible 
under v, the last equation means that W has invariant vector qj ,; so that W = Id 
is a solution for all j. Choosing this solution: (m) = 1 k=i ( i,k ^ j) for all j 

and m, means that blocked users are sent back to their original station. It can be 
interesting to note that for each j the set of solutions of <53 is convex, so that 
one can take convex combinations - with coefficient possibly depending on to - of 
different solutions. In the reversible Q case, the following rerouting is then possible: 
blocked users at end of route [ij] flip a coin, which probability of giving “head” may 
depend on the current distribution m of the other M — 1 bikes, and according to 
the result, either take route [ji], or choose new route [jk] with probability qjk- 

Other solutions can be given in particular cases, as if for example Q is uniform, 
that is, qij = (J\ — l)^ 1 for all i,j with i ^ j. Here, variants of the deterministic 
rerouting from [ij] to [ ji] can be mentioned, such as rerouting from [ij] to [jij m ], 
where to is the state of the network excluding the blocked bike, and for all j, to, 

• either i* rn is deterministic, and i i —> i* m is a one-to-one mapping on \{y}, 

• or ij m is the first step of some symmetric random walk - which distribution 
may depend on j and to - on some graph with vertex set J\ and constant 
degree. 

We will now consider any dynamics that leads to stationary distribution given 
by (15.311 . here on state space S C N M . In order to use the asymptotic results of section 
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14.21 we here again consider a sequence of networks indexed by N, and use the same 
notations as in the infinite capacity case of section 15.11 
A set of natural hypotheses is for example: 

(HO) 3 C £ N such that c^jv < C for all j £ , 

(HI) 3 yu_|_, /z_ > 0 such that /x_ < p, a ,N < P+ for all N and a £ U J™ , 

(H2) 3 A > 0 such that max v, y < A/ for all N. 

(Recall that 55 ^i,NQij,N = 1/2). Theorem l4.2l then has the following corollary. 

[ij]eJ 2 N 

Corollary 5.2. Assume that (HO), (HI) and (H2) hold and that (Mjy) satisfies 
lim Mjv =+oo and lim = 0. 

N—>oo N—too 

Then as N goes to infinity, at stationarity, the different queue-lengths at stations 
and routes get asymptotically independent, with respective approximate distributions 

• truncated-geometric with parameters Cj^jy and jNrj,N for station j, 

• Poisson with parameter TivHijLiv f or route [ij ], 

where 7jv solves equation (14.101) . Moreover, 7jv has same order of magnitude as 
M n . 


Proof. Equation (14.10|) together with the elementary relations 


°p,c 


< cp 


v c 1 o n 


V o n V o 

Z^ra=0 0 Z_m=0 0 

imply the following inequalities, using (HO) and (HI): 


< cp, for p > 0 and c £ N, 


7 n 


In 


9 < 55 ^ jvr fo],JV < M n < C 55 7 Nrj, N + 55 ^ r [ij], N < 9 

Mejfi j&Jfi Mejfi 


(C+l), 


that show that 7 jv is of the same order as M N . 

Now assumptions (1), (2) and (3) of Theorem 14.21 are satisfied, for any K > 1. 
Indeed, (1) is equivalent to (HO). As for (2), 


J N 


> 55 7iVT-[y],jV > 

M eJfi 


7 N 

2p + ' 


where 7 n > 2 p_Mn/{C + 1) goes to infinity with N. And for (3), using the 
previous lower bound on bjy, together with (H2) and inequality 7 n < 2 p + Mjy, 


b N r ^ Nr MH < 


v / 2/ct + 7iy 
2p- 


max Vi N q zj N 

imj 2 N 


^ A 0+ VMn 

- A iz~ir 


which tends to zero as N goes to infinity by assumption. 

The proof is then complete, using Theorem 14.21 □ 


Corollary 5.2 provides simple explicit approximations that may help measuring 
the performance of a real network. As an example, the total stationary rate of 
failure is given (here removing index N) by 

T = 55 = 0 ) + 55 (£[i?] • 
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Using asymptotic independence and approximate distributions as stated in Corol¬ 
lary 5.2, one gets 


E 




(7 rj) c 


= E 


(nj r An + E ^fo']7 r [ij] y^C 3 - / Y 

jeJi ^«=o!7 jJ [yjejs 2-jn—0\ l'j 1 ieji 

where we have used JT M[*j]Hy] = D, ^<7^/2 = zy/2 = /Xj7y. 


Mi 


1 + (7 r ’j) C3+1 

En=o( 7 r i)" 


If all traffic parameters are known, one can then minimize r over 7, which 
amounts to choosing the best possible M, since M and 7 are related by the one- 
to-one relation (14.101) . Since the optimal 7, and hence M, should be of order J\ (if 
rj’s are of order 1/Ji), conditions of Corollary 5.2 are satisfied. 


6. Appendix 

Proof of Theorem \2.1l 

Following the lines of the proof of the Lindeberg Central Limit theorem given in 
j2j (Theorem 27.2 p. 359), it is not difficult to show that (12.31) is satisfied, for some 
sequence (An) converging to infinity, provided that the following reinforcement of 
the Lindeberg condition is satisfied: There exists some sequence of positive real 
numbers en such that lim7v->oo £iv = 0 and 
1 AN) 

(6- 1 ) nFL IF E E - m i,w) 2l lbA—= °- 

N 3 = 1 

Now it is easily proved that the Lyapunov condition EH implies existence of 
such a sequence (en)- Indeed, assuming that (12.11) is satisfied, then the following 
inequality holds for any positive e, 


Y AN) 

-gr r E E ~' 


3,nY 1 


\Xj,N—mj } N\>ebN 


) 


AN) 


< 


J2 E (\ X TN~m j , N \ 2 + s ) 

°N j= l 


ON_ 


where S is as in (12.11) and cxn d — b N ^ 2+S) FjYi^ E (\^j,N — rn j,N\ 2+S ) , so that 

lima at = 0. It results that EH holds for en defined as En = a^ 2S \ which clearly 
satisfies lirnejv = 0. 


Proof of Theorem \2.2\ 

First note that m implies the following, apparently stronger, property 


( 6 . 2 ) 


lim sup 


2 E 


it S N~ a N. 

e b N 


- 1 


= 0, 


for some sequence (Bn) of positive real numbers converging to infinity. Indeed, 

it s "- a N ' 


assuming that (12.31) holds and setting /3 n = sup E 

that lim /3n = 0, equation (16.21) is then satisfied by any sequence (Bn) such that 


- e-‘ / 2 


so 


0 < Bn < An and e 


b 2 n /2 


<j8 


-1/2 

N 


for all N , 
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since, if these inequalities hold, then 


sup 

\t\<B N 


"E 




< e 


sup 



<13 


1/2 

N ■ 


The numbers Bn = min(Ajv, log/3jv) thus satisfy (16.21) together with limi3jv = oo. 


The proof is now standard: Using inverse Fourier transform, it results from 
E(e itSN )=J2 ¥ ’( s N = k)e itk (t e R), 


ke z 


that for any k € Z, 


P (S N = k) = ^~ [ e~ itk E(e It5iV )dt= 1 


rnb N 


27r 


27T&7V J-nb N 


-Uzn (k) Jg 


i s N~ a N 


dt. 


where we have set ZN(k ) = (k — a^/bN- 

It can be assumed without loss of generality that Bn < 7t6jv for all N > 1. 
Otherwise, replace the sequence (Bat) by (Bn A 7t&jv), which still satisfies (16.21) and 
goes to infinity, since lirnjv dv = oo by assumption (1). Then for all N and k, 


2nb N V(S N = k) - V2^e~ z »( k)/2 


' —irftjv 

/ Bn 

' 

-Bn 


e -it* N (fc) E / ) dt- 


-itz N (k)-t 2 /2 I e t 2 /2 E ( e 


-itz N (k)-t 2 /2 ^ 


— 1 dt 



[ e -itz N (k)-t 2 /2 dt 

JB N <\t\<TTb N V / 

J\t\>B N 


< 


sup 




e‘ /2 E 

4>(t)dt 


v S N ~a N 


-')L 

e~ t2 / 2 dt. 


2 / 2 dt 


l\t\>B N 


Here assumption (3) has been used to dominate the second term. Since the three 
terms of the last sum do not depend on k and converge to zero as N goes to infinity, 
the theorem is proved. 


Proof of Proposition \2.1\ 

One can assume without loss of generality that all the Poisson parameters Xj,N 
(N > 1, j € jfj) are bounded above by some fixed number, say 1 for example. 
Indeed, if this is not the case, one can replace each Poisson variable Xj n with 
A j n > 1 by the sum of [_A.j,ivJ + 1 independent Poisson variables with the same 
parameter AyjvdVwJ + l)” 1 < 1- This change does not affect Sn, &n and bN 
(only J(N) is increased). Thus the proposition will be true in the general case if it 
is proved to hold when all Ayjv are less than 1. 

Let us now show that conditions (2) and (3) of Theorem 12.21 are satisfied (note 
that (1) coincides with (ii) of the proposition). To prove that (2) holds, we use 
Theorem 12.11 and check that the Lyapunov condition is satisfied with 5 = 1. The 























26 


CHRISTINE FRICKER AND DANIELLE TIBI 


following domination of the centered third moment of any positive random variable 
X by its non centered third moment is useful: 

(6.3) E (\X - E(X)\ 3 ) < E{X 3 ). 

This results from the inequality, due to positivity of X, and hence of m = E(A'), 
\X — m \ 3 < (X + m )(X — m) 2 = X 3 — mX 2 — m 2 X + m 3 , so that E(|X — m| 3 ) < 
E(A 3 ) - mE(X 2 ) < E(X 3 ). 

Recall that the third moment of a geometric random variable with parameter p 
is given by (p + 4 p 2 + p 3 )/ (1 — p) 3 , while that of a Poisson variable with parameter 
A is A 3 + 3A 2 + A. Using the fact that all pyjv’s and Ayjv’s are less than 1, one gets 


1 V^iWl’V |3\ ^ 1 ( E {\Xj,N m j,N | 3 ) 

Wt — ^ (|x " - mj ’ Nl ) - w r [ - ZTZ - 


< 


b 3 N \i<j<J(N) 

1 E (*|iv) 

-— max -— 

bN i <j<JW 


a j,N 


J(N) 

£4 

3 =1 


N 


< 


max 


1 + + p 2 N 


bN \jejf 1 - Pj,N 


+ max (A 2 ^ + 3Xj N + 1) ] <t— ( -—-— 

jejf J b N \l-p 


where for the last step, assumption (i) is used. Then by (ii), the Lyapunov condition 
is satisfied with <5=1. 

Now to check (3) of Theorem 12.21 recall that the characteristic function of a 
geometric random variable X with parameter p is given by 


E(e itx ) = - 


1 - P 
— pe lt 

while if X is Poisson with parameter A, then 

E(e itx ) = e x{eit ~ l) . 

One easily derives that 

itx m P 


|E(e itx )| = 


\Jl — 2pcosf + p 2 
1 - P 


= 1 + 4 


(1 ~pf 


sin 2 (f/2) 


\/(l - p) 2 + 2p(l - cos t) 
in the first case, while in the second 

|E ) | = g —A(i— c °si) _ e -2Asin 2 (t/2)^ 

Then, using the convexity inequality | sin(f/2)| > |t|/7r for |f| < 7r, |E (e itSw ) | = 


J(N) 

n 

3 = 1 


Yl |E [e itXj ' N ^ | < 




n f ^ 4 


Pj,N 


7 T 2 (i - pj, N y 


exp 


£ x jx 


3£J} 


-o £ log(l+4 


je+f 


Pj,N 


7T 2 (1 - Pj,Nf 


- 2 - 


£ X J,N 




= exp 


for \t\ < 7r. 
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Due to concavity of the log, the inequality 
A t 2 Pj,N 


A P def 

< 4--— = r 


7T 2 (1 - Pj, N ) 2 (1 - P) 2 

that holds for all t £ [—7r, 7r] and all j £ J™ gives 


log ( 1 + 4 


#2 

£ Pj,N 
7T 2 (1 - pyiv) 5 


> 4 


log(l + r) t 2 Pj,N 

T 7T 2 (1 - Pj,a) 2 ' 


One obtains 


E ( e its„)| < exp _2 


log(l + r) t 


2 E 


Pj,A 


T 77" “ , (1 ^ Pj,n) 2 7T 


“ 2 r2 E 


3ZJi 


< exp —2 


log(l + r) b%f 


for |f] < 77, 


since —————- < 1. The condition (3) of Theorem 12.21 is thus satisfied with 


integrable (f> given by <j>{t) = exp ( —2 


log(l + r) t 2 


Proof of Proposition 1 2.21 

As in the preceding proof, one can assume that all Xj.N are less than one, and 
then proceed to check conditions (2) and (3) of Theorem 12.21 Condition (2) is 
obtained, here again, by proving that the Lyapunov condition holds with <5 = 1. 
This goes exactly as in Proposition ^. II except for the truncated geometric variables 
Xj, n for j £ Here instead of using (16.31) . we use the following inequality (since 

| Xj : N — mj : N | is dominated by cy at): 

E (| Xj tN - mj,jv| 3 ) < c j}N E ([X jtN - mj, N ] 2 ) = Cj,Ncr 2 N < Ca 2 N . 

This gives 


^ ,, if E (\Xj n — m,' at | 3 ) 

y E ( \Xj N - m.j N \ 3 ) < — max ——— 5 


h 3 

b N j=1 


a j,N 


b N V i<i<./(.v) 

< -— ( C + max (A 2 ^ + 3\j n + 1) J ^ -— (C + 5) 

b N V izJi J b N 


and by (ii), the last quantity goes to zero as N goes to infinity. The Lyapunov 
condition, and hence (2), is satisfied. 

As for condition (3), using Lemma [2.21 gives 


E( e itSiv)| < exp _ 2k . 


(2 E 

jejf 


J j,N 


- 2 - 


E A uV 




< exp ( —2 min(l, k ) ^ N f 


for Itl < 77. 
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